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ABSTRACT 

We present results based on an XMM-Newton observation of the prototype narrow-line Seyfert 1 
galaxy I Zw 1 performed in June 2002. The 0.6 — 12 keV spectrum is well described by a steep power-law 
(r ~ 2.3) and a weak soft blackbody component (kT ~ 200 eV) below 2 keV. The soft X-ray excess 
emission is featureless and contributes only ~ 10% to the total X-ray emission in the 0.6 — 2 keV band. 
There are tentative evidences for an iron iron Ka line and an iron K edge, both from high ionized He-like 
iron. The 0.6 — 12 keV continuum of I Zw 1 can be explained in terms of a blackbody (kT ~ 40 eV), 
thermal comptonization and its reflection off the surface of an ionized accretion disk. We note that 
characterizing the soft excess copmonent 'above a power law continuum' may not be appropriate due to 
the low energy curvature of a comptonized spectrum. X-ray emission from I Zw 1 was highly variabile 
during the observation. I Zw 1 showed a large, symmetric X-ray flare during which the X-ray luminosity 
increased by AL ~ 6 x 10 43 erg s _1 in an interval of ~ 2800 s only. The X-ray spectral variability was 
unusual during the flare. In contrast to the general observation of spectral steepening with flux from 
Seyfert galaxies, there is a clear evidence for hardening of the 0.6 — 10 keV spectrum of I Zw 1 with 
increasing flux. This hardening is due to a stronger and slightly flatter power-law component at higher 
flux without any change in the soft X-ray excess component. The large flare, the accompanying spectral 
flattening , and the steady nature of the soft excess component may be due to an additional flatter 
power-law component during the flare possibly arising as a result of ejection of electron plasma and the 
SSC process. 

Subject headings: galaxies: active - galaxies: individual: I Zw 1 - galaxies: Seyfert - X-rays: galaxies 



1. INTRODUCTION 

A significant fraction of type 1 active galactic nuclei 
(AGN) show 'soft X-ray excess emission' above a power- 
law continuum, usually identified as the steepening of the 
X-ray continuum below ~ 2 keV. This soft excess emission 
was first observed by HEAO-1 (Pravdo et al. 1981), and 
EXOSAT (Arnaud et al. 1985; Singh et al. 1985). The 
number of AGN with soft X-ray excess emission or 'ul- 
tra soft AGN' increased dramatically with the launch of 
ROSAT. Boiler, Brandt, & Fink (1996) showed that AGN 
with steepest soft X-ray spectra in the ROSAT band tend 
to lie at the lower end of the H/3 line width distribution. 
These AGN with FWHM H/3 £ 2000 km s _1 are classified 
as the narrow-line Seyfert 1 (NLS1) galaxies (Osterbrock 
& Pogge 1985), and are distinguished from the bulk of 
the Seyfert 1 galaxies ('broad-line Seyfert Is' or BLSls). 
NLSls also have strong optical Fe II emission. Detailed 
X-ray study of NLSls was carried out with ASCA which con- 
firmed the soft X-ray emission and, in addition, showed that the 
2—10 keV continuum slope too is steeper and anti-correlated 
with the FWHM of the H/3 line (Brandt, Mathur, & Elvis 1997; 
Leighly 1999b; Vaughan et al. 1999b). ASCA observations also 
showed that many NLSls have Fe Ka line arising from highly 
ionized iron, and hence the accretion disks of NLSls must be 
ionized (Dewangan 2002, and references therein; see Ballan- 
tyne, Iwasawa, & Fabian 2001 for ionized reflection model for 
NLSls). 

During the ROSAT and ASCA era, NLSls drew attention of 
most workers due to the extreme variability of their strong soft 
X-ray emission. NLSls are the most extreme X-ray variable 
objects among the radio-quiet AGN. NLSls frequently exhibit 
rapid (doubling time scales of a few hundred seconds) and/or 
large (up to a factor of 100) amplitude X-ray variability (e.g., 



Boiler, Brandt, & Fink 1996; Dewangan et al. 2001,2002). The 
excess variance for NLSls is typically an order of magnitude 
higher than that observed for BLSls. (Leighly 1999a; Turner 
et al 1999). The strong soft X-ray emission of NLSls is usu- 
ally attributed to the intrinsic emission from viscous accretion 
disks (Pounds et al. 1995). However, the variability of the soft 
X-ray emission has been found to be faster than the charac- 
teristic time scales for an accretion disk (Boiler et al. 1993). 
The fastest time scales for an accretion disk are the thermal 
time scale, t t h — 4/3q~ 1 lj~ 1 , on which thermal instabilities 
grow, and the sonic time scale, t s = r/coh, on which the sound 
waves cross the disk. In the above, a is the viscosity parameter 
(Shakura & Synyaev 1973), lo is the angular velocity and h is 
the height scale of the disk. These time scales are about a fac- 
tor of two larger than the observed variability time scale for the 
NLS1 galaxy IRAS 13224-3809 (Boiler et al. 1993). The same 
is true for many NLS1 galaxies. This apparently argues against 
the accretion disk origin for the strong soft X-ray emission and 
the nature of the soft X-ray excess emission remains uncertain. 

Chandra and XMM-Newton grating spectroscopy of NLSls 
have revealed that the soft X-ray excess emission is a 
smooth continuum component rather than a blend of emis- 
sion/absorption features (Turner et al. 2001; Collinge et al 
2001; Gondoin et al. 2002; Page et al. 2003). It is not clear 
whether the rapid variability of the smooth soft X-ray emission 
is due to the variability of the soft excess emission above and 
over the power-law or it is the underlying power-law component 
which is varying rapidly. Based on a 10 day ASCA observation 
of IRAS 13224-3809, Dewangan et al. (2002) suggested that the 
rapid variability in the soft X-ray band is likely to be due to 
the rapid variability of the power-law emission and there may 
not be intrinsically rapid variability in the soft excess emission 
above the power-law. This observation, if confirmed further, 
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will have profound effects on constraining models for X-ray 
emission of NLSls. However, this requires time resolved X- 
ray spectroscopy with short time scales which was not possible 
with the earlier X-ray observatories. XMM-Newton having the 
largest collecting area and broad bandpass is the most suitable 
for the above purpose. 

In this paper, we present results based on an XMM-Newton 
observation of I Zw 1, including time resolved spectroscopy 
on time scales of ~ 5000 s, comparable to the dynami- 
cal time scale for a Keplerian accretion disk td yn = 3.2 x 
1o3 (^) 3/2 (to^) s - 1 Zw 1 is the Prototype NLS1 galaxy. 
Its optical spectrum shows exceptionally narrow Balmer emis- 
sion lines (FWHMh/3 ~ 1240 km s _1 ; Boroson & Green 1992). 
I Zw 1 is well known for its strongest optical and UV Fe II 
emission (Boroson & Green 1992; Laor et al. 1997). Its high 
absolute luminosity (My = —23.8 for Ho = 50 km s _1 Mpc -1 , 
qo = 0) makes it the nearest quasar (z = 0.060875 ± 0.000223) 
and is also called as the 'narrow-line quasar' - a high luminosity 
version of the NLS1 galaxies. 

2. OBSERVATION AND DATA REDUCTION 

I Zw 1 was observed by XMM-Newton observatory (Jansen 
et al. 2001) on 22 June 2002 using the European Photon Imag- 
ing Camera (EPIC), the reflection grating spectrometer (RGS), 
and the optical monitor (OM) instruments. The EPIC MOS 
(Turner et al. 2001a), and PN (Striider et al. 2001) cameras 
were operated in small, and large window modes, respectively, 
to alleviate pileup of photons from the source. The EPIC cam- 
eras used the medium optical blocking filter. 

The raw events were processed and filtered using the most 
recent updated calibration database and analysis software (SAS 
v5 . 4 . 1) available in October 2003. Events in the bad pixels and 
those adjacent pixels were discarded. Only events with pattern 
0-4 (single and double) for the PN and 0-12 for the MOS were 
selected. Examination of PN light curve extracted from the full 
field but above 12 keV showed flaring particle background in 
two observation intervals of about 500 s. These intervals were 
excluded from the rest of the analysis. This resulted in 'good' 
exposure time of 18.0 ks for the PN and 18.7 ks for the MOS 
cameras. 

Source photons were extracted from the PN and MOS data 
using a circular region of radius 40" centered at the source 
position. Associated background photons were extracted from 
source free regions near the source. The total number of net 
source photons is 1.43 x 10 5 detected with the PN and 4.3 x 10 4 
detected with each of the MOS cameras. 

3. TEMPORAL ANALYSIS 

X-ray light curves of I Zw 1 were extracted from the PN 
and MOS data using a circular region of 40" centered at the 
source position and in the soft (0.4 — 2 keV), hard (2— 10 keV), 
and full (0.4 — 10 keV) bands. The soft and hard bands were 
chosen to separate approximately the two spectral components 
- a power-law and a soft excess component generally observed 
from NLS1 galaxies (see e.g., Leighly 1999b). The time bins 
with sizes of 200 s are at least 50% exposed in all the light 
curves. Background light curves were extracted from source 
free regions with the same bin sizes, exposure requirements, 
and in the same energy bands as for the source light curves. 
The background light curves were subtracted from the respec- 
tive source light curves after appropriate scaling to compensate 
for the different areas of the extraction regions. The three light 
curves extracted from the PN and MOS data, were combined 
to improve the signal-to-noise ratio. The observed counts cor- 
respond to a mean observed flux of 1.9 x 10~ n erg cm~ 2 s _1 , 
and luminosity of 1.6 x 10 44 erg s _1 in the 0.4 — 10 keV band. 

Figure 1 shows the combined PN and MOS light curves of 
I Zw 1 in the full, soft, and hard bands. Also shown in Fig 1 



is the hardness ratio defined as the ratio of count rates in the 
2 — 10 keV and 0.4 — 2 keV bands. It is evident that X-ray emis- 
sion from I Zw 1 varied strongly during the XMM-Newton ob- 
servation. The light curves show trough-to-peak variations by a 
factor of ~ 1.6 in the full band, ~ 1.5 in the soft band, and ~ 2 
in the hard band. The most remarkable variable event is a large 
X-ray flare which started at ~ 10 4 s after the beginning of the 
observation and lasted for about 5000 s. During this flare the 
full band count rate increased from 10.16 ± 0.14 counts s" 1 to 
14.43 ±0.15 counts s _1 in just 2800 s. The count rates quoted 
above were calculated as the mean of at least three data points. 

The soft and hard band light curves show similar variability 
pattern and appear to be correlated. To quantify this correla- 
tion and to search for any possible time lag, we have calculated 
cross-correlation function (CCF) of the hard band light curve 
with respect to the soft band light curve. For this purpose we 
utilized PN and MOS light curves with time bins of 16 s created 
in a similar way as described above. The CCF was calculated 
using the FTOOLS program crosscor. The CCF is plotted in 
Figure 3 with time bins of 256 s, which shows strong correla- 
tion between the soft and hard band light curves without any 
significant delay. We also examined CCFs with smaller time 
bins and found no time lag between the two light curves. 

Despite the strong correlation between the soft and hard 
band light curves, the X-ray flare has different amplitude of 
variability in the soft and hard bands. The flare is more in- 
tense in the hard band (an increase in count rate by a factor 
of ~ 1.6) than in the soft band where the count rate increased 
by a factor of ~ 1.4. This is reflected in an increased hardness 
ratio during the flare (see Fig. 1). It appears that the hardness 
ratio generally follows the variations in the count rate in any 
band. This is obvious in Figure 2 where the hardness ratio is 
plotted against the count rate in the full band. Figure 2 shows 
that the X-ray emission from I Zw 1 hardens with X-ray in- 
tensity. The trend in the hardness ratio with count rate tells 
us that the total counts in the 2—10 keV band increases more 
than that in the 0.4 — 2 keV and nothing about the 2 — 10 keV 
or 0.4 — 2 keV spectral change. 

4. TIME-AVERAGED SPECTRAL ANALYSIS 

Photon energy spectra of I Zw 1 and associated background 
spectra were accumulated from the EPIC PN and MOS data 
using the same source and background extraction region de- 
scribed above. The pulse invariant channels were grouped such 
that each spectral bin contained at least 80 and 30 counts in the 
case of PN and MOS spectra, respectively. The EPIC responses 
were generated using the SAS tasks rmfgen and arfgen. 

The first order RGS1 and RGS2 spectra of I Zw 1 were ex- 
tracted using source regions in dispersion/spatial coordinates 
that contained 90% of the point source events. The correspond- 
ing background spectra were generated using the entire regions 
of CCD arrays but excluding 95% of the point spread function 
along the spatial direction. RGS response matrices were cal- 
culated using the SAS task rgsrmf gen. Both the RGS spectra 
were binned to 5 channels, a further grouping (a minimum of 5 
counts per spectral channel) was applied to avoid zero counts 
in any channel. 

All the spectral fits were performed with the XSPEC pack- 
age (version 11.2; Arnaud 1996) and using the x 2 -statistic for 
the EPIC spectra and C-statistic for the RGS spectra. Unless 
otherwise specified, the quoted errors on the best-fit model pa- 
rameters were calculated for the 90% confidence level for one 
interesting parameter i.e., A\ 2 or AC = 2.71. 

4.1. EPIC spectra 

At first the MOS and PN spectra were fitted separately to 
check for the possible uncertainties due to cross-calibration 
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problems. We found generally good agreement between MOS2 
and PN cameras (Ar ~ 0.03) in the 2.5 - 10 keV bands, how- 
ever, MOS1 results in a flatter (Ar ~ 0.15) spectrum. These 
discrepancies are consistent with that found by Molendi & Sem- 
bay (2003). In the 0.6 - 2.5 keV band, MOS1 spectrum is flat- 
ter (Ar ~ 0.13) than PN spectrum, while MOS2 data agrees 
well with the PN data (Ar ~ 0.04). Therefore we present the 
spectral results obtained by fitting the same model jointly to 
the PN and MOS2 data while allowing the relative normal- 
izations to vary. However, PN and MOS data do not agree 
below ~ 0.6 keV. Similar discrepancies at low energies have 
been found for PG 1211+143 (Pounds et al. 2003) and Akn 
564 (Vignali et al. 2003). Therefore, we did not use EPIC data 
below 0.6 keV for out spectral analysis. 

In order to characterize the continuum shape of I Zw 1, first 
we fitted a redshifted and absorbed power-law model jointly to 
the PN and MOS2 data. Por this fit we chose the energy bands 
2.5 — 6 keV and 8 — 10 keV bands in the observer's frame. These 
bands are relatively free of any spectral features in the AGN 
spectra. This model resulted in a good fit (\ 2 = 304.7 for 292 
degrees of freedom (dof)). The best-fit power-law photon index 
is very steep (r~x ~ 2.28). A comparison of the observed data 
with the extrapolated best-fit power-law to lower energies is 
shown in Figure 4 in terms of deviations x = {N° bs — N* nod )/a i , 
where Nf is the observed counts in energy channel i, cr, is 
the standard error on iV° bs , and A 7 l mod is the best-fit model 
counts in channel i. Fig. 4 shows a weak soft X-ray excess 
emission above the power-law at energies below ~ 2 keV, a 
weak line-like emission at ~ 6.7 keV, and an edge-like feature 
at ~ 8 keV. The soft X-ray excess emission is a general feature 
of NLS1 galaxies. To model this soft component we added a 
blackbody component to the power-law model and carried out 
the fitting in the 0.6 — 10 keV band. This model provided a 
good description of the spectra (x 2 = 851 for 785 dof). The 
best-fit absorption column density is Nn = 8.31 2 ,' 2 , x 10 20 cm~ 2 , 
which is significantly higher than the Galactic column density 
of N§ al = 5.0 x 10 20 cm" 2 (Dickey & Lockman 1990) along the 
line of sight to I Zw 1. This suggests an additional absorber lo- 
cal to the source. Therefore, we fixed the column density to the 
Galactic value and added an additional photoelectric absorp- 
tion model at the source redshift to our two component model. 
This model provided a good fit (x 2 = 846.3 for 785 dof) re- 
sulted in an excess absorption of AA^h = 5.llJ'g x 10 20 cm~ 2 . 
As also noted earlier, there are signatures of possible presence 
of an iron Ka line and K-edge. Addition of a redshifted absorp- 
tion edge to the model consisting of a blackbody and a power- 
law modified by Galactic and intrinsic absorption improves the 
fit (Ax 2 = —9.3 for 2 additional parameters). Adding a Gaus- 
sian line further improves the fit (x 2 = 822.1 for 780 dof). The 
iron K-edge and the iron Ka line are detected at significance 
levels of 98.6% and 99.7%, respectively, according to an F-test. 
Figure 5 shows the observed PN and MOS2 data, the best-fit 
model, and the deviations x (see above) of the data from the 
best-fit model. 

4.2. RGS spectra 

Due to its excellent spectral resolution, the RGS instrument 
allows a more detailed examination of the soft X-ray spectrum. 
The RGS spectra of I Zw 1 are shown in Figure 6. The soft 
X-ray spectrum of I Zw 1 does not show strong emission or 
absorption features. It is a smooth continuum except for a 
few weak features. Therefore, we fitted an absorbed power-law 
model jointly to the RGS1 and RGS2 spectra in the 0.35—2 keV 
band. The relative normalizations were allowed to vary. The 
model resulted in C-statistic = 857.5 for 805 dof with signif- 
icant intrinsic (excess over the Galactic column) absorption 
of AA^h = Q-St^'.o x 10 20 cm~ 2 and a steep power law with 
r = 2. 411^20 • The excess absorption is similar to that de- 



rived from the PN and MOS data. However, the soft X-ray 
power- law appears to be slightly steeper (Ar ~ 0.1) than the 
hard X-ray (2.5 — 10 keV) power-law derived from the PN and 
MOS data. The best-fit model and ratio of the observed data 
and model are shown in Fig. 6. There are no strong sharp fea- 
tures. The only obvious features are in the 0.5 — 0.6 keV band, 
probably associated with the detector response of the oxygen 
K-edge, and emission line-like feature at a rest-frame energy 
of ~ 1.9 keV. Adding a narrow Gaussian line at ~ 1.94 keV 
results in AC = —7.3, and an equivalent width of ~ 20 eV. 
The rest-frame line centroid energy is not well constrained but 
appears to be higher than the centroid energy of Si Ka line 
(1.74 keV for neutral Si). 

4.3. The soft X-ray excess emission 

The blackbody and power-law model described the data ad- 
equately, but the temperature of the blackbody (kT ~ 230 eV) 
is unrealistically high for a standard accretion disk (see Section 
8.4). Appropriate accretion disk models in place of a simple 
blackbody result in similar temperatures. The multicolor disk 
blackbody model diskbb provided equally good fit (x 2 = 821.3 
for 780 dof) and resulted in kT = 309±^j eV, while the more 
realistic disk model diskpn (Gierlinski et al. 1999) provided 
equally good fit and resulted in kT = 300^49 eV. 

4.4. Ionized reflection model 

The presence of an iron Ka line at ~ 6.7 keV and an iron K- 
edge at ~ 8.5 keV both suggest that these features arise from 
highly ionized iron (Fe XXIII - Fe XXV). The strength (EW 
~ 130 eV) and width of the iron line (FWHM ~ 27000 km s _1 ) 
suggests that the line is produced in the accretion disk. These 
observations are supportive of Compton reflection emission 
from a highly ionized accretion disk. The reflection emission 
from ionized material can be strong below ~ 2 keV and con- 
tribute significantly to the soft X-ray excess emission. With 
this in mind we adopt an X-ray continuum model that in- 
cludes a power-law and Compton reflection from ionized ac- 
cretion disk. We used the Compton reflection model pexriv of 
Magdziarz & Zdziarski (1995) which calculates Compton re- 
flection of illuminated X-ray power-law continuum by an ion- 
ized accretion disk. The pexriv model does not calculate the 
iron Ka line, therefore we also included a Gaussian line for 
the iron Ka line. As before, we used the photoelectric ab- 
sorption due to Galactic and intrinsic columns. This model 
resulted in a good fit (% 2 = 831.1 for 781 dof). The best-fit 
parameters were F = 2.521°;^, R = 1.6l° g (where R = 1 
corresponds to a solid angle of 2n subtended by the reflector), 
£ = 773^773* erg cm s _1 . Slightly higher value of R may sug- 
gest an additional spectral component for the soft X-ray excess 
emission. Addition of a blackbody component further improves 
the fit. The model consisting of a blackbody, power- law, re- 
flection and a Gaussian, and modified by Galactic as well as 
intrinsic absorption (model 2) improves the fit (x 2 = 820.8 
for 779 dof). The best- fit parameters were F = 2.41±gg| , 
R = 0. 651^7, £ = 1683t?^ 4 erg cm s" 1 , and parameters of 
the Gaussian line were E = 6.72±°H keV, a = 225±225 eV > 
and EW = 89 (see Table 1). The best-fit model and the un- 
folded spectrum are shown in Figure 6. Due to the contribution 
of the ionized reflection in the soft X-ray band, the strength of 
the blackbody component is weaker in model 2 compared to 
that derived in model 1. 

4.5. Ionized reflection and thermal comptonization 

Inverse Compton scattering of soft photons by hot elec- 
trons, possibly located in a corona above the disk, can pro- 
duce approximate power-law continuum. The primary contin- 
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uum emission of AGN is a comptonized spectrum not an ex- 
act power law. The departure at low energies of a thermally 
comptonized spectrum from a power-law continuum may have 
profound effects on the shape of the soft X-ray excess emission 
of NLS1 galaxies. In model 2, the soft X-ray excess compo- 
nent is modeled by ionized reflection and a weak blackbody 
component. Since the thermal comptonized spectrum is an ap- 
proximate power-law, we further investigate the significance of 
the weak blackbody component after replacing the power-law 
component in model 2 by a comptonization model. We used 
the thermal comptonization model comptt (Titarchuk 1994). 
The free parameters of the model are the temperature of the 
seed photons (kT s ), assumed to follow a Wien law, the optical 
depth of the scattering region (r), the electron plasma tem- 
perature (kT e ), and a parameter defining the geometry (a disk 
or a sphere). In the following, we consider the disk geometry. 
We fixed the parameters of the Gaussian line and the illumi- 
nating spectrum in the pexriv model to be the same as that 
of model 2 but varied the reflection fraction R. This model 
provided a good fit (x 2 = 815.2/783). The best- fit blackbody 
temperature is now 38^2 e V- The best-fit relative reflection 
is 0.8 ± 0.1. The best-fit parameters of the comptt compo- 
nent are seed photon temperature kT a = 155^20 e V, electron 
temperature kT e = 12.2tt% 8 keV, and r = 1.3t°; 2 . 

4.6. Two component thermal Comptonization model 

A corona with two electron populations with distinct tem- 
peratures can give rise to both a soft X-ray excess and a power- 
law. It is likely that the surface of an accretion disk is ionized, 
the outermost layers may contain thermal electrons in addi- 
tion to the hotter electron population above the disk in the 
corona. Under such conditions, the cool disk photons will be 
first comptonized by the cooler of the electron population at 
the surface of the disk, giving rise to the soft excess emission 
below ~ 2 keV. This soft excess emission may be the seed 
photons for the comptonization by the hotter electrons in the 
corona (Page et al. 2003). We investigated if such a process 
is a viable model for the X-ray spectrum of I Zw 1. First, we 
tried to fit the power-law component with the comptt model. 
A combination of a blackbody, comptt, a Gaussian line, and 
an edge model provided a good fit (\ 2 = 825.0 for 779 dof). 
No intrinsic absorption is required in this model. Thus the 
power-law component is well described by comptonization of 
a soft thermal spectrum with temperature kT = 134^Jf eV 
by hot electrons with kT e — 24.4+g'?, and the soft X-ray ex- 
cess component is described by a blackbody component with 
kTtb = 292^jy eV. The contribution of this soft excess compo- 
nent to the 0.6 — 2 keV band unabsorbed emission is ~ 11%. As 
the temperature of the soft excess component is much higher 
than that expected from an accretion disk. This component 
is unlikely to be the intrinsic disk emission. Replacing the 
blackbody compoent by an additional comptonization compo- 
nent the fit slightly (x 2 = 820.7 for 777 dof). The best-fit 
parameters are comptt 1: soft photon kT — 60^43, electron 
temperature kT = 16. 51 28 ! 2 and r = 1.71^2 i comptt 2: soft 
photon kT = 19lt 2 g eV, electron temperature kT = 29.8t.tH 
and t = 0.29to 5 16 . In this model, soft thermal photon from 
the disk with kT = 60±g are comptonized by hot thermal 
(kT = 16. 51 28 ! 2 ) electrons in the corona. This process gives 
rise to the soft X-ray excess emission. The second comptoniza- 
tion model produces the power-law emission. 

5. TIME-RESOLVED SPECTROSCOPY 

To examine the spectral variations of I Zw 1, we extracted 4 
time-selected spectra across the 20 ks XMM-Newton observa- 
tion. We chose the sampling time scale of 5000 s, which sepa- 
rates the flaring interval nicely. Thus we extracted four spectra 



- spectrum 1 (0 - 5000 s), spectrum 2(5001 - 10000 s), spec- 
trum 3 (10001-15000), and spectrum 4 (15001-20000 s). The 
resulting "on-source" average exposure time was ~ 5 ks per 
EPIC spectrum. We set the ancillary response, and response 
matrix files to be those of the mean spectrum. Since there 
is no variation in the background level across the observation, 
we used the same background spectra extracted across the full 
observation for the mean spectra. Again, the energy channels 
were appropriately grouped to achieve a good signal-to-noise 
(a minimum of 30 and 20 counts per energy channel for the PN 
and MOS, respectively). As before, we carried out the spectral 
fitting of the time-selected PN and MOS2 data jointly. 

To illustrate the flux and spectral variability of I Zw 1, we 
compared each of the time-selected spectra with the mean spec- 
trum. Figure 9 shows deviations of the time-selected spectra 
from the best-fit model 1 for the mean spectrum. These plots 
were constructed by using the best-fit model 1 for the mean 
spectrum but without actually fitting the time-selected spectra. 
Strong variation is present only for the flux. To quantify these 
variations we performed spectral fitting of the time selected 
spectra. We used the model 1 (a combination of a blackbody 
+ power law + Gaussian line multiplied by Galactic, intrin- 
sic absorption and an iron K-edge models) and fitted the 4 
time-selected EPIC spectra. It was not possible to constrain 
both the Gaussian line and iron K-edge for individual spectra. 
Therefore, we fixed the line energy and width of the Gaussian 
line and the parameters of the edge model to be the same as 
the best-fit values derived for the mean spectrum. Model 1 
provided good fit to all the the four time-selected spectra of 
I Zw 1. The best-fit parameters are listed in Table 2. The 
results of the time-resolved spectroscopy are plotted in Fig- 
ure 10, which shows that all the parameters except for the 
power-law normalization changed significantly across the 20 ks 
XMM-Newton observation. The decrease in the ratio of the flux 
of the blackbody and power-law components in the 0.6 — 2 keV 
band during the large flare appears to be due to increase in the 
power-law flux alone. 

To investigate further, we fit the model 1 to the time-selected 
PN spectra jointly. We excluded the MOS2 data from this fit to 
avoid large number of model parameters. As before, the width 
and the centroid energy of the iron Ka line, and the parame- 
ters of the edge model were fixed to the best-fit values derived 
from the mean spectrum. In this fit, except for the normaliza- 
tion of the power-law component, all the parameters for a data 
group were tied to the corresponding parameters in other data 
groups. The four normalizations of the power-law components 
for the four time-selected were varied independently. This type 
of parameterization allowed us to check if the spectral vari- 
ability of I Zw 1 is solely due to variations in the power-law 
normalization. The joint fit to the time-selected spectra re- 
sulted in a good fit (x 2 = 1889.1 for 1903 dof). Varying the 
photon indices independently improves the fit (Ax = 49.1 for 
3 less dof). There is a marginal evidence for flattening of the 
power-law during the large flare. The best-fit photon indices 
are 2.36±g;g|, 2.33±g;g|, 2.24±g;g£, and 2.29t^l for spectrum 
1, 2, 3, and 4, respectively. Varying blackbody temperatures in- 
dependently does not improve the fit (x 2 = 1834 for 1897 dof). 
The same is true for the blackbody normalization (x = 1835.4 
for 1897 dof). 

6. FLUX-SELECTED SPECTROSCOPY 

The soft and hard band X-ray emission are correlated with- 
out any significant delay. Based on our time-resolved spec- 
troscopy in the preceding section, this correlation appears to 
be due to the variability of the power-law component alone 
without any significant changes in the soft X-ray excess com- 
ponent. In order to further explore and to improve the de- 
pendence of the spectral parameters, we have carried spectral 
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anaiysis at iow and high flux ievels. We extracted two aver- 
aged spectra covering the large flare between the elapsed time 
interval of ~ 10 4 — 1.5 x 10 4 s (high flux or flare state) and 
excluding the large flare (low flux or quiescent state). The 
mean PN count rates are 6.7 ±0.03 count s -1 (high flux state) 
and 5.2 ± 0.02 count s -1 (low flux state) in the 0.6 - 10 keV 
band. We analyzed both the PN spectra jointly. We used the 
model 1 (absorbed and redshifted blackbody and power-law 
modified by an absorption edge at ~ 8.5 keV and iron Ka line 
at ~ 6.7 keV. We fixed the parameters of the iron line and 
edge to be the same as that derived for the mean spectrum. 
In order to investigate which model parameter or parameters 
differ most significantly between the two flux states, we tied all 
the model parameters for the high flux state to that of low flux 
state and untie the parameters step-by-step. 

The best-fit model with the same model parameters for the 
two flux states is not a good fit (x 2 = 2782.8 for 1281 dof). 
Varying the power-law normalizations independently resulted 
m an acceptable fit (x 2 = 1347.8 for 1280 dof). Untying 
the photon indices further improved the fit (\ 2 = 1302.3 for 
1279 dof), with significantly different best-fit photon indices of 
T = 2.24± J ; )4 (high flux state) and 2.34+° ° 4 5 (low flux state). 
We also varied the blackbody normalization independently, the 
fit did not improve further (\ 2 = 1300.0 for 1278 dof). Thus 
we conclude that the blackbody normalization did not vary be- 
tween the two flux states and tied the parameter again in the 
following analysis. We calculated the confidence contours for 
the photon indices and the power-law normalizations in the low 
and high flux states. Figure 11 shows these contours at 68%, 
90%, and 99% confidence levels. 

The brightening and flattening of the power-law during the 
X-ray flare suggests that the flare could have resulted due 
to an additional power-law component. To test such a sce- 
nario, first we tied the photon index and normalization for 
the two flux states and introduced a second power-law com- 
ponent for the high flux state. This model resulted in a good 
fit (x 2 = 1305.2 for 1279 dof). The best-fit photon index and 
the normalization of the second power-law are F = 2.011q;q7 
and npL = 1.55t.o Qg x 10~ 3 photons cm -2 s _1 , respectively. 



The mass is then 



M = 1.464 x 10' 



The first power-law is steeper and stronger (r 
ii PL = 5.36t{5;26 x 10 -3 photons cm -2 s -1 ). 



9 ,,+0.08. 



7. BLACK HOLE MASS FOR I ZW 1 

Before we discuss the results below, we estimate the mass of 
the black hole for I Zw 1. Calculations of the black hole mass 
rely on the assumption that the dynamics of the optical broad- 
line region gas is dominated by the gravitational potential of 
the central super-massive black hole. To calculate the black 
hole mass for I Zw 1, we used the results of Kaspi et al. (2000), 
based on a reverberation study of 17 quasars. Kaspi et al. 
(2000) have determined an empirical relationship between the 
size of the broad-line region (Rblr) and the monochromatic 
luminosity (\L\) as follows. 



Rblr = (32.9) 



AL A (5100 A) 



10 44 erg s- 



lt - days (1) 



We use the flux taken from a linear interpolation between pho- 
tometric data points obtained from the Third Reference Cat- 
alog of Bright Galaxies (RC3; de Vaucouleurs et al. 1991) 
at the rest wavelength of 5100 A (A/ A (5100 AA) = 6.1 x 
10 -11 erg cm -2 s -1 ). The mass of the black hole is given by 
Mbh = rv 2 G~ 1 . To determine v, the velocity, we correct 
vfwhm of the H/3 emission line by a factor of V3/2 to ac- 
count for velocities in three dimensions We used vfwhm = 
1240 km s -1 for the H/3 line from Boroson & Green (1992). 



Using the cosmology go = 0.5 and Hq = 70km s _1 Mpc _1 , this 
gives us a value for the central black hole mass of 2.4 x 1O 7 M . 

8. DISCUSSION 

This paper presents an analysis of ~ 20 ks XMM-Newton ob- 
servation of the prototype narrow- line Seyfert 1 galaxy I Zw 1. 
The time-averaged X-ray spectrum of I Zw 1 shows four spec- 
tral components - a power-law continuum, a soft X-ray ex- 
cess component below ~ 2 keV, weak intrinsic absorption 
(Nh ~ 2 x 10 20 cm -2 ), and spectral features due to iron (the 
Ka line and K-edge). Above 2.5 keV, the power-law compo- 
nent has a slope of F ~ 2.3 which is steeper than that of most 
BFS1 galaxies (r ~ 1.9) but is similar to other NLS1 galaxies 
e.g, Ton S180 (r = 2.26l°'?2; Vaughan et al. 2002), Mrk 335 
(r = 2.29±0.02; Gondoin et al.2002); Akn 564 (F ~ 2.50-2.55; 
Vignali ct al. 2003) . The X-ray continuum slightly steepens at 
lower energies, resembling to the soft X-ray excess emission 
generally observed from NLS1 galaxies. The strength of this 
soft X-ray emission is weaker in I Zw 1 than that in many 
NFS Is. The contribution of the soft X-ray excess emission to 
the 0.6 — 2 keV band X-ray emission is only 10% in I Zw 1. 
ASCA observation of I Zw 1 in 1995 also showed intrinsic ab- 
soption (Nh ~ 3 x 10 20 cm -2 ), steep power-law continuum 
(rx ~ 2.4), and an iron Ka line at ~ 7 keV but did not detect 
any soft excess component (Leighly 1999). I Zw 1 was fainter 
by a factor of ~ 2 during the ASCA observation. 

8.1. Rapid variability and the efficiency limit 

True to the nature of NFS1 galaxies, I Zw 1 showed rapid 
and large amplitude X-ray variability during the XMM-Newton 
observation. X-ray emission from I Zw 1 increased by 40% in 
2800 s only, corresponding to a change in the 0.6 — 10 keV 
luminosity of AL = 6 x 10 43 erg s _1 . This is a remarkable 
variability for a radio-quiet quasar as the time scale is a factor 
of ~ 6 shorter than the dynamical time scale for a 10 7 M© 
Schwarzschild black hole. Such variability is rare among radio- 
quiet quasars. Other radio-quiet quasars to show such rapid 
and large amplitude variability are PHF 1092 (Brandt et al. 
1999), RX J1334.9+3759 (Dewangan et al. 2002). 

The major consequence of rapid and large amplitude vari- 
ability is that the source has to be emit very efficiently. It is 
possible to calculate the efficiency of conversion of rest mass 
into energy using simple arguments originally due to Fabian 
(1979). If the X-ray emitting region is an accreting spheri- 
cal cloud of ionized hydrogen, then there is a limit on the ra- 
diative efficiency, n, given by n > 4.8 x 10 -43 AZ//Ai, where 
AL is the change in luminosity in the rest-frame time in- 
terval At. First, we convert the change in the combined 
PN+MOS count rate into change in the luminosity by using 
a constant conversion factor of (1.5 ± 0.1) x 10 43 erg count -1 , 
derived from the time averaged PN+MOS count rate and 
the intrinsic luminosity. The change in the count rate of 
3.99 ± 0.14 count s -1 over 2800 ± 565.7 s (2639.3 ± 533.2 s 
in the rest frame) corresponds to a rate of change of lumi- 
nosity, AL/ At = (2.27 ± 0.49) x 10 40 erg s -2 . This gives a 
lower limit to the efficiency, n > 0.011 ± 0.002. This efficiency 
can be compared with the maximum efficiency possible for a 
Schwarzschild black hole (n ~ 0.06; Shapiro & Teukolsky 1983) 
and a Kerr black hole (n ~ 0.3; Thorne 1974). Thus the rapid 
flare observed from I Zw 1 is consistent with that expected for 
a Schwarzschild black hole. 

Our time resolved spectroscopy has revealed that the soft 
X-ray excess emission, presumably arising from an accretion 
disk, did not vary appreciably during the giant flare seen in 
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Fig. 1. The rapid and large amplitude variability is consistent 
with a flaring of the power-law component alone. The primary 
power-law emission of radio-quiet AGN is thought to be the 
result of comptonization of soft photons from the disk by hot, 
relativistic electrons in a region, the geometry of which is not 
clear at present. The possibilities are - (i) an extended corona 
above a cold accretion disk (e.g, Haardt & Maraschi 1991), (ii) 
non-thermal electrons at the base of a jet, and (Hi) thermal 
electrons in the innermost regions of an accretion disk. In the 
first case, rapid and large amplitude variability can be produced 
by magnetic flares arising from the surface of the disk. In the 
second case, ejection of a relativistic plasma blob as a jet would 
produce the giant flare with power-law spectral shape. In the 
third scenario, it is difficult to produce giant and rapid flares. 
It is difficult to distinguish between the first two possibilities 
in the absence of wide band X-ray spectrum up to several hun- 
dred keV. However, a Kerr black hole is not required in any 
of the two possible scenarios. Magnetic flares can be produced 
on very short time scales as the reconnection time scale can be 
very short (Merloni & Fabian 2001). 

8.2. Flux correlated spectral variability 

There is a clear evidence for a change in the hardness ratio 
with total count rate (see Fig. 2) . The full band X-ray spectrum 
of I Zw 1 hardens with increasing flux. A detailed time-selected 
spectral analysis revealed that hardening of the full band spec- 
trum as suggested by the hardness ratio is due to the increase 
in the power-law normalization and a steady soft excess com- 
ponent. Our flux-selected spectroscopy has revealed that the 
power-law component became stronger and slightly flatter in 
the flaring or high flux state. Thus the change in the hardness 
ratio is not only due to the variations in the 2—10 keV band flux 
but also due to hardening of the primary power-law. The spec- 
tral hardening of I Zw 1 with increasing flux is in sharp contrast 
with that observed from many Seyfert galaxies, which gener- 
ally show spectral steepening with increasing flux (Turner 1987; 
Singh et al. 1991; Turner et al. 1990; McHardy, Papadakis & 
Uttley 1998; Lee et al. 2000; Chiang et al. 2000; Done et 
al. 2000; Zdziarski & Grandi 2001; Petrucci et al. 2001; Shih, 
Iwasawa & Fabian 2002; Vaughan and Edelson 2001; Nandra 
2001; Romano et al. 2002; Dewangan et al. 2002; Gallo et 
al. 2003; Lamer et al. 2003). Georgantopoulos & Papadakis 
(2001) found similar behavior in a number of Seyfert 2 galaxies. 
To explain such spectral variability, McHardy et al.(1998) and 
Shih et al. (2002) proposed a two component model, consisting 
of a soft component with constant spectral shape but variable 
flux and a hard component with constant flux, so that the soft 
component dominates the spectrum at high fluxes and overall 
spectrum becomes steeper. Using a long XMM-Newton ob- 
servation, Fabian & Vaughan (2003) showed that the spectral 
variability of MCG-6-30-15, on time scales of 10 ks, is charac- 
terized by a soft varying component which is a power-law with 
constant slope and a hard constant component produced by 
very strong reflection from a relativistic disk. In contrast, in 
the case of I Zw 1, it is the soft component with constant flux 
and the hard component with variable flux resulting in spectral 
hardening. 

8.3. Possible origin of spectral hardening 

The observed 2—10 keV spectral steepening with flux has 
been attributed to the change in the shape of the comptonized 
spectrum possible due to Compton cooling by increased seed 
or soft photon flux from the disk (Haardt, Maraschi, & Ghis- 
ellini 1997) or due to increasing covering fraction of the corona 
resulting from correlated magnetic flare or reconnection events 
('thunder-cloud' model of Merloni & Fabian 2001). Obviously 
both the scenarios are not applicable to the spectral behavior 



of I Zw 1. The soft excess emission did not increase with flux 
during the flare, implying that there was no change in either 
the seed photons from the disk or reprocessed emission by the 
disk. The soft excess emission attributed to the intrinsic disk 
emission and reprocessed emission simply did not appear to 
respond the large X-ray flare. This is possible if the flare emis- 
sion was beamed away from the accretion disk. The X-ray flare 
and accompanied spectral hardening of I Zw 1 is suggestive of 
a generation of a new population of hot or relativistic electrons 
without an increase in the covering fraction of the corona. This 
new population must be short-lived and beamed away from the 
disk. This is possible if the new population is ejected from the 
base of a jet and may be part of a discrete jet. The X-ray 
flaring then could be due to synchrotron self-Compton (SSC) 
process as in the case of radio-loud AGN. It is interesting to 
note the similarities of the differences between the X-ray spec- 
tra of radio-loud and radio-quiet AGN and that of I Zw 1 in 
the low flux and high flux states. 

8.4. The soft X-ray excess emission 

The high resolution soft X-ray spectrum of I Zw 1, obtained 
with the RGS, is featureless and smooth. We did not detect 
any emission or absorption features below 1.5 keV, the region 
of the spectrum usually affected by the warm absorbers. Thus, 
the RGS spectrum of I Zw 1 supports the earlier observations 
with Chandra LETG and XMM-Newton RGS that the soft X- 
ray excess emission of NLS1 galaxies is featureless and smooth 
(Turner et al. 2001; Collinge et al 2001; Gondoin et al. 2002; 
Page et al. 2003). Although the soft X-ray excess emission 
above a power-law continuum is well described by a blackbody 
component, its temperature is too high for the intrinsic emis- 
sion from an accretion disk. The temperature of a standard 
thin accretion disk is 

/ \ -3/4 

T(r) ~ 6.3 x i0 5 m 1/4 M 8 - 1/4 ( -LJ K (3) 

where rs is the Schwarzschild radius, and Ms is the black hole 
mass in units of 10 s M© (Peterson 1997). The temperature at 
the last stable orbit (r — 3rs) for a Schwarzschild black hole is 
~ 35 eV for a black hole mass of 2.4 x 10 7 M© accreting at the 
Eddington rate. This is an upper limit to the disk temperature, 
because m < 1, and outer regions of the disk will have lower 
temperatures. 

A more realistic model involving ionized reflection in ad- 
dition to the blackbody and power-law components describes 
the 0.6 — 12 keV spectrum well. The reflection component fits 
the iron edge and contributes significantly to the soft excess 
emission thereby reducing the strength of the blackbody com- 
ponent. However, the blackbody temperature is still too high 
for an accretion disk. 

The primary X-ray continuum of Seyfert type AGN is 
thought to arise from thermal Comptonization of soft photons 
from the disk by the hot electrons in the corona (Haardt & 
Maraschi 1991). Such a comptonized spectrum deviates from 
the power-law form at low energies, in addition to the high 
energy cut-off, due to the Maxwellian distribution of hot elec- 
trons (see e.g., Titarchuk 1994). Therefore, power-law is not 
an appropriate continuum to characterize the soft X-ray ex- 
cess component. In the case of I Zw 1, the soft X-ray excess 
emission above a 'thermally comptonized continuum' is well de- 
scribed by ionized reflection and blackbody. Most importantly 
the temperature of the blackbody (kT ~ 40 eV) is consistent 
with that expected for the inner regions of an accretion disk 
around a ~ 10 7 M© black hole. Thus the soft X-ray emission 
of the proto-type NLS1 galaxy I Zw 1 is not unusual in terms 
of strength, temperature or variability time scale. 
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8.5. Iron features 

We detected an iron Ka line at ~ 6.7 keV and an iron K- 
edge at ~ 8.5 keV. Both the features arise from He-like iron 
and suggest that the accretion disk of I Zw 1 is ionized. Ev- 
idence for ionized accretion disks have been found in many 
NLS1 galaxies e.g., Ton S180 (Vaughan et al. 2002; Comastri 
et al. 1998) Mrk 335 (Gondoin et al. 2002), Akn 564 (Vaughan 
et al. 1999) and are consistent with the general picture that 
accretion disks of AGN with steep X-ray spectrum are ion- 
ized (Dewangan 2002). The ionization stage of photo- ionized 
material is determined by the ionization parameter, £ = -^i, 
where L is the luminosity of ionizing photons, n is the density 
of the material, and R is the distance of the material from the 
ionizing source (Kallman & McCray 1982). For standard thin 
disks, Matt et al. (1993) showed that the ionization parameter 
strongly depends on the accretion rate relative to the Edding- 
ton rate, £ oc {-^ — ) 3 . In order to ionize He-like iron, £ should 

be at least ~ 500 or j±— > 0.2 (Matt et al. 1993). This result 
is consistent with the best-fit value of the ionization parame- 
ter derived from ionized reflection model. These observations 
have important implications to the nature of the soft X-ray 
excess emission as reflection from ionized material contributes 
substantially in the soft X-ray band, as revealed in the case of 
I Zw 1. 

9. CONCLUSIONS 

We presented timing and spectral characteristics of the 
proto-type NLS1 galaxy I Zw 1 based on an XMM-Newton 
observation. The main results are as follows. 

1. The 0.3 — 12 keV spectrum of I Zw 1, obtained with 
XMM-Newton, consists of four intrinsic spectral compo- 



nents, namely a steep (Fx ~ 2.3) primary continuum, a 
soft X-ray excess component described by a blackbody 
(kT ~ 40 eV) and ionized reflection from the disk, iron 
features - Ka line and K-edge both from He-like iron 
arising due to reflection of the primary continuum, and 
an intrinsic photoelectric absorption component. 

2. I Zw 1 showed a large X-ray flare with changes in the 
luminosity AL ~ 6 x 10 43 erg s _1 on a time scale of 
~ 2800 s. The radiative efficiency, n = 0.01 ± 0.002, 
inferred from the flare is consistent with X-ray emis- 
sion outside the last stable orbit around a Schwarzschild 
black hole. 

3. X-ray spectrum of I Zw 1 hardened with flux during the 
observation. This behavior is uncommon among Seyfert 
galaxies. The hardening of the 0.6 — 10 keV is not only 
due to increase in the hard (2 — 10 keV) band flux but 
also due to hardening of the power-law component at 
higher fluxes. The soft excess emission did not respond 
to the changes in the power-law continuum. The spec- 
tral variability of I ZW 1 is difficult to understand in 
the framework of disk-corona models but appears to be 
related with ejection of relativistic plasma at the base of 
a jet. Simultaneous X-ray and radio observations would 
be required to understand such spectral variability dur- 
ing large X-ray flares. 

This work is based on observations obtained with XMM- 
Newton, an ESA science mission with instruments and contri- 
butions directly funded by ESA Member States and the USA 
(NASA). REG acknowledges NASA award NAG5-9902 in sup- 
port of his Mission Scientist position on XMM-Newton. 
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Fig. 1 — Combined PN and MOS light curves of I Zw 1 in the 0.4 - 10 keV, 0.4 - 2 keV, and the 2-10 keV bands. The bottom panel 
shows hardness ratio defined as the ratio of count rates in the 2 — 10 keV and 0.4 — 2 keV bands. 



0.26 



0.24 



s. 

w 
w 



0.22 



T3 

I 0.2 



0.18 



i 1 1 1 1 1 1 1 1 1 1 1 r 

I Zw 1 





_l I I L 



_l I I I I L 



10 12 14 

Count rate (count s -1 ) 

Fig. 2. — Hardness ratio Vs the 0.4 — 10 keV count rate showing spectral hardening with increasing flux. 
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Fig. 3. — Cross-correlation function for the soft (0.4 — 2 keV) and hard (2 — 10 keV) band EPIC light curves of I Zw 1 with 256 s bins. The 
two bands are strongly correlated without a delay. 
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Fig. 4. — A comparison of the observed data in the 0.6 — 10 keV band with the best-fit power-law (r ~ 2.28) model derived using the data 
in the 2.5 — 6 keV and 7—10 keV bands, and extrapolated at lower energies. The deviations are defined as \ = (N° hB — N* nod )/<Ti, where 
N° ha is the observed counts in energy channel i, cr; is the standard error on N° bs , and 7V™ od is the best-fit model counts in channel i. A 
weak soft X-ray excess emission below ~ 2 keV and a weak iron a line is seen at ~ 6.7 keV 
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Fig. 5. — PN and MOS2 spectral data and the best-fit model consisting of an absorbed and redshifted blackbody (kT ~ 235 eV) and 
power-law (r ~ 2.3), a Gaussian emission line at ~ 6.7 keV, and an absorption edge at ~ 8.5 keV. 
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Fig. 8. — Unfolded PN and MOS2 spectra and the best-fit model consisting of a (1) blackbody (kT ~ 40 eV), (2) reflection from an ionized 
accretion disk, (3) a Gaussian line, and (4) thermal Comptonization model (comptt). 
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Fig. 9. — Deviations of the individual time-selected PN (open circles) and MOS2 (filled circles) spectra from the best-fit model to the mean 
spectrum. The time-selected spectra are (a) spectrum 1, (b) spectrum (2), (c) spectrum 3, and (d) spectrum 4. 
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Fig. 10. — Time series of the best-fit spectral parameters derived from the time-resolved spectral fittings. The time bin size is 5000 s. 
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Fig. 11. — Confidence contours of the power-law photon index and normalization for the low and high flux state spectra. The contours are 
drawn at the confidence levels of 68%, 90%, and 99%, and the best- fit values are marked by a plus sign. 
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Table 1 

Best-fit spectral model parameters derived from the joint spectral fitting of the PN and MOS2 spectra of 



I Zw 1. 








Spectral model a 


Comp. 


Parameter 


BB+PL+G+Edge 


BB+PL+G+Refl BB+G+CompTT+Rcfl 






(model 1) 


(model 2) (model 3) 



Excess Absorption 
Blackbody 



Power law 



CompTT 



Gaussian 

Edge 

Pexriv 

Total 



AN H (xlCH" cm" 2 
kT ( eV) 

n b b b 

f bb (0.6 - 2 keV) f 

r 

n p i c 

f P l{0.6-2 kcV) f 
/ p; (2-10keV) f 
kT see( i eV 
kT p i keV 

T 

n c d 

/c(0.6 - 2 keV) f 
fc(2 - 10 keV) f - 

a°(cV) 
ft- s 

J line 

EW (eV) 

E e dge 
T 

R 

£ ( erg cm s -1 ) 
/ o6a (0.6-10 keV) f 
f mtr 0.6 - 10 keV) f 
xLiJdof 



234.5±^7 

q0.5 

0.11 

9 97+0.05 
z - z ' -0.05 
+0.35 



5.56 



0.39 

0.92 
0.88 



u - ' z -0.20 
2fi7 +307 
ZD ' -149 

i o +0 - 7 

llU -0.4 

130 

°- oo -0.84 

0.14lg;» 



1.6 
1.9 
822/780 



8.5Z22 
199^9 

1 c+0.9 

0.08 

9 41+O.O6 

z - 41 -0.09 

R O+0.46 

1.0 

0.85 



0.75 



+0.22 
0.20 
+ 112 
-225 
n y+0.6 
U -'-0.5 

89 



6.72 
225 



0.65+°f 
1683l?°66 

1.6 
2.0 
820.8/779 



z - z -1.6 

38l« 



2278 



+2.5 xlO 5 
-2268 



0.007 



1^^-20 
i 9 9+44.8 
iz - z -6.3 

1 6+ 31 
1 -°-0.3 
i q+0.2 

0.9 

6.72 
225 
0.7 



o.83i8:» 

1683 
1.6 
1.9 
815.2/783 



a BB - blackbody; PL - power law; G - Gaussian; Edge - absorption edge; Refl - ionized reflection model (pexriv); CompTT - Thermal 
comptonization model. The spectral models were modified by the Galactic (N^ al = 5 X 10 20 cm -2 ) and intrinsic absorption (ANh). 
^Blackbody normalization in units of 10~ 5 X 10 39 erg s -1 /(d/10 kpc) 2 , where d is the distance. 
c Power-law normalization in units of lO - 3 photons cm -2 s~ 1 keV - 1 at 1 keV. 
d CompTT normalization in units of 10 — 3 photons cm -2 s _1 keV -1 at 1 keV. 
e Error on this parameter was calculated after fixing the edge parameters at their best-fit values, 
■fin units of 10 — 11 erg cm -2 s -1 
9 Iron line flux in units of 10 
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Table 2 

Results of time resolved spectral fitting of the PN and MOS2 spectra of I Zw 1 









BB+G+PL a 




Comp. 


Parameter 


Spectrum 1 


Spectrum 2 


Spectrum 3 


Spectrum 4 


Excess absorption 


A7V H (xlO 20 cm" 2 




6-2+3^ 




3-413.3 


Blackbody 


kT ( cV) 


244.5t^-? 


201.1^-? 


298.3+^-3 


275.8_ 63 2 




n b b d 


9 7 +0.9' 
z -'-1.0 


o 1+3.1 
°- i -2.2 


2.0tif 


9 i+l.O 
z - i -0.8 




M0.6 - 2 keV) b 


0.13 


0.13 


0.11 


0.13 


Power law 


r 


9 9O+0.09 
z - zo -0.08 


9 01 +0.08 
z - ol -0.08 


9 9 r+0.08 
z - zo -0.09 


9 99+O.O9 
z - zz -0.10 




n p i c 


4 - 8 -0.6 


c 9+0.5 
°- z -0.6 


6.8185 


4 6 +0 - 7 

4 - D -0.8 




/ pi (0.6-2keV)s 
f p i(2 - 10 keV)s 


0.79 


0.85 


1.12 


0.76 




0.71 


0.74 


1.06 


0.74 


Gaussian 


Er d 

inline 


6.71(/) 


6-71(/) 


6.71(f) 


6.71(/) 




<T d (eV) 


266.5(/) 


266.5(/) 


266.5(/) 


266.5(/) 




f h 

J line 


7 7+8.6 
'■'-7.2 




10.8+^ 


5.3±g-J 




EW (eV) 


117 


231 


107 ' 


76 


Edge 


E e dge 


8.64(/) 


8.64(/) 


8.64(/) 


8.64(/) 




T 


0.19(/) 


0.19(/) 


0.19(/) 


0.19(/) 


Total 


/ o6s (0.6-10 keV)s 


1.46 


1.50 


2.04 


1.45 




/ mtr 0.6 - 10 keV)s 


1.64 


1.73 


2.31 


1.62 




X 2 rmJdof 


588.2/612 


556.2/600 


726.6/731 


709.3/648 



a A combination of a blackbody (BB), power-law (PL), Gaussian line (G), and an absorption edge (Edge) modified by the Galactic absorption 
as well as intrinsic or excess absorption (ANn). 

^Blackbody normalization in units of 10~ 5 X 10 39 erg s -1 /(d/10 kpc) 2 , where d is the distance. 
c Power-law normalization in units of 10 -3 photons cm -2 s -1 kcV -1 . 

d The line energy and the Gaussian a were fixed to their best- fit values derived from the mean spectrum of I Zw 1. 

f Fixed parameters 

9 In units of 10 -11 erg cm -2 s — 1 

h Iron line flux in units of 10 — 6 



